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achieve more versatile, convenient stability criteria for a wide class of 
finite-difference approximations to initial boundary value problems associated 
with the hyperbolic system u t = Au^ + Bu + f in the quarter-plane 
x >_ 0, t >_ 0. With these criteria, stability is easily established for a 
large number of examples, thus incorporating and generalizing many of the 
cases studied in recent literature. 
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0. Introduction. 


In this paper we sharpen and extend the results of [4] in order to achieve more ver- 
satile, convenient, sufficient stability criterion for a large class of approximations to the 
initial-boundary value problem associated with the hyperbolic system 
u f = Au, + Bu + f in the quarter plane x > 0, t > 0. Our difference approxima- 
tion consists of a general difference scheme -- explicit or implicit, dissipative or not, 
two-level or multi-level -- and boundary conditions of quite a wider type than discussed 
in [4]. 

As in [4], we restrict attention to the case where the outflow components of the 
principal part of the boundary conditions are translatory; i.e., determined at all boun- 
dary points by the same coefficients. Such boundary conditions are commonly used in 
practice; and in particular, when the boundary consists of a single point, the boundary 
conditions are translatory by definition. 

Throughout the paper we assume that the basic scheme is stable for the pure Cau- 
chy problem, and that the other assumptions which guarantee the validity of the 
Gustafsson-Kreiss-Sundstrom stability theory in [5], hold for our case. With this, we 
raise the question of stability for our approximation in the sense of Definition 3.3 in [5] . 

Our stability analysis begins in Section 2 where we show (Theorem 2.1) that our 
entire approximation is stable if and only if the scalar outflow components of its princi- 
pal part are stable. Thus, our global stability question is reduced to that of a scalar, 
homogeneous, outflow problem which, as in [4j, is the main subject of this paper. 

We state our stability criteria for the reduced problem in Theorems 3.1 and 3.2 of 
section 3. These criteria depend both on the basic scheme and the boundary condi- 



tions, but very little on the intricate interaction between the two. It follows that our 
criteria provide in many cases a convenient, easily-checkable alternative to the well- 
known Gustafsson-Kreiss-Sundstrom criterion in [5]. 

We proceed in Section 3 to use our stability criteria in Theorems 3.1 and 3.2 
together with Lemmas 3.1 and 3.2 in order to establish all our previous examples in [4] 
as well as new ones. This includes a host of dissipative and nondissipative examples 
that incorporate and generalize many of the cases studied in recent literature; e.g. [1]- 
|10|, [12]-[15|. 

As in [4] we point out that there is no difficulty in extending our stability criteria 
to two-boundary problems, since if the corresponding left and right quarter-plane prob- 
lems are stable then, by Theorem 5.4 of [5], the original problem is stable as well. 

We also remark that there are no essential obstacles in extending our results to 
initial-boundary value problems with variable coefficients. 
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1. The Differential Problem and the Difference Approximation. 

Consider the first-order hyperbolic system of partial differential equations 

(1.1a) du.(x,t)/dt = Adu(x,t)/dx + Bu{x,t) = f {x,t), x > 0, t > 0, 

where u(ar,l) = (uW(x, /),•••, u^ n \x,t)Y is the unknown vector (prime denoting the 
transpose), T{x,t) = {fl l \x,t),--‘, jl n \x,t))' a given n-vector, and A and B fixed 
n X n matrices such that A is symmetric and nonsingular. Without restriction we 
may assume that the system is given in characteristic variables, namely A is diagonal 
of the form 


A 1 0 
0 A n 

where A 1 and A n are of orders i X € and (n -€) X (n -€), respectively. 
The solution of (1.1a) is uniquely determined if we prescribe initial values 

(1.1b) u(ar,0) = u(ar), x > 0 , 

and boundary conditions 



, A 1 > 0, A n < 0 , 

j 


(1.1c) .u n (0,f) = Su ! (0,<) + g(t), t> 0, 

where S is a fixed (n-€) X € coupling matrix, g (t) a given (n-t)- vector, and 

(1.3) u 1 = (u^V"j «^)'» u 11 = (u^ +, V") 

a partition of u into its outflow and inflow components, respectively, corresponding 
to the partition of A in (1.2). 
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Introducing a mesh size Ax > 0, At > 0, such that X = At/ Ax = constant, 
and using the notation v^t) = v(i/Ax,t), we approximate (1.1a) by a general, con- 
sistent, two-sided, solvable basic scheme of the form 


+ At) = £ Q ff y„{t -a At) + A/b„(<) 

< 7=0 


i/ = r, r + !,••• 


(1.4) 


Qff ^2 Aj Ev„ v^+j, <r — 1, ,s, 

;'=-r 

where the n X n coefficient matrices ^ are polynomials in XA and A<5, and the 
n-vectors !>„(/) depend on f (x,t) and its derivatives. 


The difference equations in (1.4a) have a unique solution vjl + At) if we pro- 
vide initial values 


(1.5) 


v„(/iAt) = v^/xAf), /i=0, •••,«, u = 0,1,2, — , 


and specify, at each time level t — fiAt, n = s, s + l,-*- , boundary values 
v„(< + At), v = 0, “, r - 1. These boundary values will be determined by boundary 
conditions of the form 


T { -1 v v {t + At) = £ v„(< -<rA/) + A/d„(Z) , v = 0,- -, r -1 , 

< 7=0 

(1.6a) 

JlT» = E CM E\ r = -1, -, 9> 

y=o 


where the n X n matrices depend on .A, A IB and S; and the n-vectors 

d u {t) are functions of f (x,Z), g(f), and their derivatives. 

We shall assume that the leading coefficients are nonsingular, thus assuring 
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that the boundary conditions (1.6a) can be solved for the required boundary values 
v u (t + At), v = r - 1, — , 0, in terms of neighboring values of v^. 

We shall further assume that the matrices Cj£ depend weakly on 5, in that B 
introduces a mere 0(At) perturbation in these matrices. This assumption holds for 
all practical boundary conditions where the elements of are polynomials in the 
entries of AtB. 

Finally we assume that in accordance with the partition of A in (1.2), the Cj%) 
can be written as 


(1.6b) 

a. 

ii 

'c)i cyw ' 

C“IW ennw 


where 




(1.6c) 

the C)i are independent of v , 

(1.6d) 

the Cjl are diagonal when 

5 = 0, 

(1.6e) 

the C}”W 

= 0 when B = 

= 0, 

and 




(1.61) 

emu) = o 

for v > r - j 

when 5 = 0 


Let us now set 




f \ 

V 1 

w l > 


dl 

(1.7a) 

V «/ = 

V 11 
V V 

/ 


d" 
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according to the partition of u in (1.3); so that the boundary conditions (1.6a,b) split 
into two groups: 

r|.[vi(f + At) + + At) 

(1.7b) = £ [rj'vj(/ -trAl) + -<rA() I + A(d[(() , 

< 7 = 0 


T" = E Cjl&' K nM - E <$<’>&, cr = -1, „ = 0,-,r - 1; 

y =0 ;=0 


and 


T[\’Myl(l + At) + r!J ,,M v“(( + At) 

(1.7c) = £ [rJ'MvJll -aAt) + rJ“Mv['(f -<rAf)| + Afd"(f) , 

<7=0 

7»4) = j c?"MF, o=I, II, ff = -l,--,q, i/ = 0,-,r-<, 

;=o 

which we refer to as the outflow and inflow boundary conditions, respectively. With 
this formulation it is a simple matter to verify that the boundary conditions in our pre- 
vious papers [3,4] constitute a special case of the present ones. Hence, the argument 
concluding Section 1 of [3] implies that boundary conditions of the form (1.6) can be 
constructed to any degree of accuracy. 

It should be pointed out that the outflow boundary conditions in (1.7b) are quite 
general despite the apparent restrictions in (1.6c)-(1.6e). Indeed, (1.6c) is not much of 
a restriction since in practice the outflow boundary conditions are often translatory, 
i.e., determined at all boundary. points by the same coefficients. In particular, if the 
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numerical boundary consists of a single point, then the boundary conditions are trans- 
latory by definition, so (1.6c) holds automatically. The restrictions in (1.6d)-(1.6e) pose 
no great difficulties either, since they are satisfied by all reasonable boundary condi- 
tions, where for B = 0 the CyJ usually reduce to polynomials in the diagonal block 
A\ and the Cj vanish. 

We realize that in view of the restriction in (1.6f), the inflow boundary conditions 
(1.7c) are not quite as general as the outflow ones. If, however, the boundary consists 
of a single point then such conditions can be achieved in a trivial manner, simply by 
duplicating the analytic condition (1.1c) which gives 

yS(t + At) = Sv l 0 (t + Af) + g (t + AO • 
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2. The Reduced Approximation. 


The difference approximation is completely defined now by (1.4)-(1.6), and we wish 
to apply to it the stability theory of Gustafsson, Kreiss and Sundstrom in [5]. For this 
purpose we assume from now on that the basic scheme (1.4a) is stable for the pure 
Cauchy problem (-00 < u < oo) and that the other assumptions in [5] are satisfied as 
well. With this, the Gustafsson-Kreiss-Sundstrom theory holds for our case* and we 
ask whether approximation (1.4)-(1.6) is stable in the sense of Definition 3.3 of [5]. 

In Theorem 2.1 which concludes this section, we shall reduce the above stability 
question to that of a scalar outflow approximation with homogeneous boundary condi- 
tions. In order to obtain this reduction we begin, as in (3], by recalling Lemma 10.3 of 
[5] which provides a necessary and sufficient determinental stability criterion stated 
entirely in terms of the homogeneous part of the difference approximation. This cri- 
terion immediately implies that for stability purposes we may consider (1.4)-(1.6) with 

(2.i) Mf)=d„(i)=o. 

Moreover, since B introduces a mere 0(At) perturbation in the matrix 
coefficients A Jcr and Cj^, then by Theorem 4.3 of [5], our approximation is stable if 
and only if it is stable for 5=0. Setting B = 0, we get (1.6d)-(1.6f), which together 
with (2.1) proves that our stability problem is equivalent to that of the principal part 
of (1.4)-( 1.6), given by a basic scheme of the form 


‘The applicability of [5] to our approximation is further discussed in the first paragraph 
of Section 2 in [3|. 
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(2.2a) 


Q- 1 v i/(* + Af) = E Qc V v( l ~ ^ = r, r + 1, •, 

(T«= 0 

Q<r= it A j<r E ’’ -1 

j=-r 


with initial values 


(2.2b) v„(/iA<) = v„(/zAf), fi = 0,-,s, v = 0,1,2,— , 

and outflow and inflow boundary conditions described, respectively, by 


ri{ vj,(< + A/) = £ r”v*(f-crA0, 1 / = (),•••, r - 1 , 

(T=0 

ri 1 = £ Cj* Ef , <r = -l ,•••,», 
y-o 

(2.2c) 

nonsingular , 


Cyj diagonal and independent of v , 


and 

rLV (,/) vj,(f + At) +T?/ I,, ' ) v“(« + Af) 

= t [7? IM »1(< -<rA<) + ryx-tyjll -<rA()|, 

•7=0 

(2. 2d) 

r ll« M = £ cM-'E', o-I,II, o = -1,-, q , 
y=o 

Cj|5f) nonsingular , z/ = (),•••, r — 1; and = 0 , z/ > r - j . 


= 9 - 



This yields: 


LEMMA 1.2 (compare 3, Lemma 2.1). Approximation ( 1.4)-( 1.6) is stable if and 
only if its principal part in (2.2) is stable. 

Since now B = 0 and b„(/) = 0 then evidently, the basic scheme (2.2a) is con- 
sistent with the principal part of (1.1a), i.e., with the homogeneous system: 


(2.3) 3u(x,f)/3f =■ Adn(x,t)/dx . 

Furthermore, B = 0 implies that the coefficients Aj a of the basic scheme (2.2a) 
are now polynomials only in \A. Thus, the A ]a are diagonal matrices; and putting 



4 0 

o 4 


in accordance with the partition of A in (1.2), we may split the basic scheme (2.2a) 
and the initial values (2.2b) into two parts: 


vl(t + At) = E -ffA<), v = r, r + 1, •• , 

( 7=0 

(2.4a) 

Ql = £ 44 <7=-l,-,< , 


(2.4b) vJ,(/iAf) = vJ,(/zAt), p, =0,-'S, v = 0,1,2,-** ; 

and 
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(2.5a) 


<?V?(I + AO = £ <?“v»(< -<rAf), !/ = r, r + 1,- , 

( 7=0 

<?? = £ a'}„e’, <t = -i, •••,«, 

/=-r 

(2.5b) vJ(aiA 0 = vJ(/iA0, A*=0,---,s, v = 0,1,2, 

We therefore view approximation (2.2) as made of outflow and inflow parts given by 
(2.4) (2.2c) and (2.5)(2.2d), respectively. And clearly, (2.2) is stable if and only if both 
parts are. 

We observe that the outflow approximation (2.4)(2.2c) is self-contained, so we can 
solve it to obtain the outflow values 

v[(f + Ai), v — 0,1,2,-" . 

Having found these values we realize, however, that the inflow boundary conditions 
(2. 2d) can now be solved for the required inflow boundary values 

(2.6) v”(f + A/), i/ = 0,-,r- 1, 

independently of the interior inflow values vj, 1 , v > r. It thus follows that the inflow 
boundary values in (2.6) can be specified in advance, merely by means of outflow 
values. Hence, the stability question for approximation (2.5)(2.2d) is equivalent to that 
of the inflow basic scheme (2.5a) with arbitrary inhomogeneous boundary values. Now, 
the above-mentioned Lemma 10.3 of [5] implies again that without affecting stability we 
may replace these arbitrary inhomogeneous values by homeogeneous ones: 
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(2.7) v?(f + Af)=0, v = 0,”,r - 1 . 

This gives us the new, self-contained inflow approximation (2.5)(2.7), whose stability 
together with that of the outflow part in (2.4)(2.2c) is equivalent to the overall stability 
of approximation (2.2). 

Since the A ]a and CyJ are diagonal, we proceed to write 

Aja = diag(a /<7 ), C) l a = diag(c Jt7 ) , 
and further split (2.4)(2.2c) into € scalar components, each of the form 


(2.8a) 


Q. \ v v{t + A<) = v = r, r + 1,— , 

< 7=0 


Qa= f; Oja®, <* = 3, 

y— r 


(2.8b) 


v„{fiAt) = v^fiAt), p = 0,‘ -,s , v =0,1,2, " , 


T_ ,v v {t + At) = £ T a v u (t -aAt), v = 0,-,r - 1 , 


( 7=0 


(2.8c) 


T„= E c y<7 ° c 0( _i) 7^ 0 , 

y=o 


where (2.8a) is consistent with a corresponding outflow component of (2.3): 


(2.9) 


du(x,t)/dt = adu(x,t)/dx, a = constant > 0 . 
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Similarly, we split the inflow approximation (2.5)(2.7) into n - ( scalar com- 


ponents: 


(2.10a) 


Q_iV„{t + At) = £ Q a v v {t -trAt), v = r, r + 1, — , 

( 7=0 


Qc= £ a io& . = 

y«=-r 


(2.10b) 


^(/zAt) = ^(//At), fi=0,-,s, v = 0,1,2, ••• , 


(2.10c) v„(t + At) = 0, v = 0,“-,r - 1 , 

where now, (2.10a) is consistent with an inflow component of (2.3): 

du(x,t)/dt = adu(x,t)/dx, a < 0 . 

Since approximation (2.1) is stable if and only if (2.4)(2.2c) and (2.5)(2.7) are, and 
since the latter are stable if and only if their scalar components are, we immediately 
obtain: 

LEMMA 2.2 [3, Lemma 2.2j. Approximation (2.2) is stable if and only if the 
approximations in (2.8) and (2.10) are stable for all eigenvalues a > 0 of A 1 , and 
a < 0 of A u , respectively. 

With this lemma in mind we proved: 
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LEMMA 2.3 [3, Lemma 2.3]. The inflow scalar approximation (2.10) is uncondi- 
tionally stable for every eigenvalue a < 0 of A u . 

This lemma -- due to Kreiss [7] in the case where the basic scheme is dissipative, 
explicit and two-level — combined with the previous two, finally yield the main result 
of this section: 

THEOREM 2.1 (compare [4, Theorem 1.1], [3, Theorem 2.1]). Approximation 
(1.4)-(1.6) is stable if and only if the reduced outflow scalar approximation in (2.8) is 
stable for every eigenvalue a > 0 of A 1 . That is, approximation (1.4)-(1.6) is stable if 
and only if the scalar outflow components of its principal part are all stable. 
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3. Statement of Main Results and Examples. 


Theorem 2.1 implies that from now on we may restrict our stability study to the 
scalar outflow approximation in (2.8). Hence, we open this section by stating the fol- 
lowing five assumptions which hold throughout this paper, and guarantee the validity 
of the Gustafsson-Kreiss-Sundstrom theory in [5j for this approximation. 

ASSUMPTION 3.1 ((4, Assumption 1.1); (5, Assumption 3.1]). Approximation (2.8) 
is boundedly solvable; i.e., there exists a constant K > 0 such that for each 
y € € 2 (Az) there is a unique solution w € f 2 (Ax) to the equations 

Q-i w u = V» v= r, r + 1, — , 

T^w v = y v , v = (),•••, r - 1 , 

with |MI<KMI; where Q_\ and T _ j are defined in (2.8a, c), and € 2 (Ax) is the 
space of all grid functions w = with ||u;|| 2 = Ax ££1^ \w„\ 2 < co. 

ASSUMPTION 3.2 ((4, Assumption 1.2); (5, Assumption 5.1]). The basic scheme 
(2.8a) is stable for the pure Cauchy problem -oo < u < oo. That is, defining the 
basic characteristic function 

(3.1) P{z,k) = £ aj(z)K J , 

y=-r 

where 

(3.2) a/z) = a /H) - £ z ^ a j<7 , j = -r,--,p , 

< 7=0 
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then: 


(i) The basic scheme fulfills the von Neumann condition ; i.e., the roots z(ac) of 
the basic characteristic equation 

P{z,k) — 0 

satisfy 

|z(/c)| < 1 /or a N K |#c| = 1 . 

(ii) If |k| = 1, and if z(/c) is a root of P(z,k ) with 1 2 (fc) | =1, then z(k) is 
a simple root of P(z,k). 

ASSUMPTION 3.3 (compare (4, Assumption 1.3], and [5, Assumption 5.4 together 
with Definition 10.1]). The basic scheme (2.8a) belongs to one of the following classes: 

(i) Dissipative schemes; i.e., schemes for which the roots z(/c) of P{z,k) satisfy 

|z(k)| < 1 for all k with |/c| = 1, k 1 . 

(ii) Almost-dissipative schemes, where the roots of P{z,k) satisfy 

]z(ac)| < 1 for all but a finite number of k with |k| =1, 

(iii) Unitary schemes (also known as strictly nondissipative schemes) where the 
roots of P{z,k) satisfy 

| z(k) | = 1 for all k with |k| — 1 . 

Obviously, if the basic scheme belongs to any of these classes, then by continuity 
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it fulfills the von Neumann condition in Assumption 3.2. 


ASSUMPTION 3.4 ([4, Assumption 1.4]; [5, Assumption 5.5]). 

a P ( 2 ) 7^ 0 f° r atl M > 1 ■ 

ASSUMPTION 3.5 [4, Assumption 1.5]. 

m 

E I c ;'( z ) I > 0 f° r al1 M > 1 > 

y=o 

where in analogy with (3.2), 

(3.3) cj(z) = c y( _ t) - £ z^~ l c ja , j 

< 7=0 

The above five assumptions are satisfied by all reasonable approximations of form 

(2.8). 

In order to state our new stability criteria we define, in complete analogy with 
(3.1), the boundary characteristic function 

R( 2 >k) = E c j( z ) Kl 
/*= o 

where the Cy(z) are given in (3.3). Putting 

U (z,k) = |P( 2 ,k)( + | /?(ar,Ac) | , 

we shall prove in Section 4, the following modification of Theorem 2.1 in [4]: 
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THEOREM 3.1 (1st Main Theorem). Approximation (2.8) is stable if 


f 2 ( 2 ,K) > 0 for all 

(3.4) 

{|z I = | K | = 1, (z,k) 5 ^ (1,1)1 U {|z | > 1, 0 < |k| < 1} . 

This result is an improvement of Theorem 2.1 in (4] in that here, the (z,k) 
domain on which Q(z,k) must not vanish is a proper subset of the corresponding 
domain in (4j. 

As in [4] we shall find it convenient to divide this (z,k) domain into three disjoint 
parts, and restate Theorem 3.1 as follows: 

THEOREM 3.1' (Theorem 3.1 restated). Approximation (2.8) is stable if 

(3.5a) f2(z,/c) > 0 for all \z\ = |k| = 1, #c ^ 1 , 

(3.5b) Q{z,k = 1) > 0 for all \z\ = 1, z 1 , 

(3.5c) Q(z,k) > 0 for all \z\ > 1, 0 < |/c| < 1 . 

The advantage of this setting over that of Theorem 3.1 will be soon clarified by 
Lemma 3.1 in which we provide helpful, sufficient conditions for each of the three ine- 
qualities in (3.5) to hold. Before stating this lemma we need, however, the following 
definitions: 

DEFINITION 3.1. The boundary conditions (2.8c) fulfill the von Neumann condi- 
tion if the roots z(/c) of the boundary characteristic function R(z,k) satisfy 

|z(k)| < 1 for all |/cj = 1 . 
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DEFINITION 3.2. The boundary conditions (2.8c) are dissipative if the roots of 
R(z,k) satisfy 


|z(/c)| < 1 for all |/c| = 1 , k 1 . 

Clearly these definitions are analogous to those made for the basic scheme in 
Assumptions 3.2 and 3.3. Further, if the boundary conditions are dissipative then they 
obviously fulfill the von Neumann condition. 

With the above definitions we can now state: 


LEMMA 3.1 (1st Main Lemma). 

(i) Inequality (3.5a) holds if either the basic scheme (2.8a) or the boundary condi- 
tions (2.8c) are dissipative. 

(ii) Inequality (3.5b) holds if any of the following is satisfied: 

(a) The basic scheme is two-level. 

(b) The basic scheme is three level and 


(3.6) n {z = - 1 , k = 1 ) > 0 . 

(c) The boundary conditions are two-level, and at least zero-order accurate as 
an approximation to equation (2.9). 

(d) The boundary conditions are three-level, at least zero-order accurate, and 
(3.6) is satisfied. 

(iii) Inequality (3.5c) holds if the boundary conditions fulfill the von Neumann con- 
dition, and are either explicit or satisfy 
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Cj(_^k 3 7^ 0 for all 0 < |/c| < 1 . 

jctsO 

This lemma follows immediately from Theorem 2.2 and Lemma 2.1 in [4). Part (iii) 
of the lemma is associated with important observations on solvability by Osher [11]. 

With Theorems 3.1, 3.1' and Lemma 3.1 at hand, it is now a trivial matter to rees- 
tablish Examples 2.1-2.10 in [4], precisely as it was carried out using Theorems 2.1, 2.1', 
2.2 and Lemma 2.1 there. Thus we conveniently achieve the following five examples: 

EXAMPLES 3.1 [4, Examples 2.1 and 2.2). Stability holds for any basic scheme of 
form (2.8a), with boundary conditions generated by either the explicit right-sided, dis- 
sipative Euler scheme: 

v u (t + At) = v v {t) + Xa[vn(*) - v„(01> 0 < Xa < 1, v = 0,-, r - 1 , 

or by its implicit analogue: 

v„(t + At) = v„(f) + \a[v v+l (t + At) - v v {t + At)], \a > 0, i/ = 0,-, r - 1 . 

EXAMPLE 3.2 [4, Examples 2.3 and 2.8). We have stability for any two-level basic 
scheme of form (2.8a), with boundary conditions determined by either horizontal extra- 
polation of order k - 1: 

»„(< + A() = £ (*) (-l)' +l tv + ,(f + Ai), v = 0,- r - 1 , 
y-i 1 

or by the right-sided, dissipative, three-level Euler scheme: 
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V v {t + AO = v u (t - At) + 2Xa[ Vfl (f - AO - v v {t - AO) , 


0 < \a < v = 0, •••, r - 1 . 

EXAMPLE 3.3 [4, Examples 2.4 and 2.5). Stability holds for any dissipative basic 
scheme of form (2.8a), with unitary* boundary conditions generated by either oblique 
extrapolation of order k - 1: 

(3.7) »„(( + At) = £ (*) (-l)> +! »„+,(( -(; -1)A(|, i/=0,-,r-l, 

/= i J 

or by the second-order accurate Box-scheme: 


v v {t + AO + v v+l (t + AO - X«k+i(< + AO - + AO] 

(3.8) 

= v„{t) + ^+,(0 + MVh( 0 - v„(0]> v = (),•••, r - 1 . 

EXAMPLE 3.4 (4, Examples 2.6 and 2.7). We have stability for any basic scheme 
of form (2.8a) whose characteristic function satisfies 


(3.9) P{z = -1, k = -1) t* 0 , 

with boundary conditions determined by the right-sided three-level, almost-dissipative, 
weighted Euler scheme: 


*The boundary conditions are unitary if the roots z(k) of R(z,k) satisfy 
1 2 (/c) | =1 for all |/c| = 1 . 
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v„{t + At) = v u {t - At) + Xa[2»„ +1 (0 - v^t + At) - v v {t - Af)j , 

(3.10) 

0 <\a <1, v = 0,*", r - 1 . 

We realize that if the basic scheme is dissipative then by Assumption 3.3(i), 

(3.11) P(z,k) 0, |*| > 1, \k\ = 1, k 7^ 1 ; 

so (3.9) holds automatically. As noted already in Example 2.7 of (4], however, certain 
well-known nondissipative schemes satisfy (3.9) as well. This includes implicit 
unconditionally-stable cases such as the unitary Crank-Nicolson scheme: 

v v {t + At) [v u+l {t + At) - v vA (t + Af)] 

(3.12) 

= v u (t) + ~ K+i(0 - Vl(01. v = r, r + 1, • , 

and the almost-dissipative, backward Euler scheme: 

(3.13) v„(t + At) - -y- K+jt* + Af) - v„_ t (t + A/)) = v v (t), u = r, r + 1,— . 

EXAMPLE 3.5 [4, Examples 2.9 and 2.10]. The Crank-Nicolson scheme in (3.12) 
with 0 < \a < 1, and the backward Euler scheme in (3.13), are stable with oblique 
extrapolation at the boundary: 

(3.14) »„(f + A!) = £ (*) - U - 1)A(] . 

/= i J 

While Examples 3.1-3.5 incorporate and extend many of the cases discussed in 
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recent literature (e.g., [l]-[8], [10], [12]-[15]), other interesting examples are not covered 
by Theorems 3.1 and 3.1'. For instance, Gustafsson et al. showed in Theorem 6.2 of [5] 
that the Leap-Frog scheme 


(3.15) 


v u {t + Af) = v(t -At) + XalvnfO - v uA (t)) , 


0<Xa<l, v = 1,2,3,"- , 

with the oblique boundary extrapolation in (3.14), provides a stable approximation. In 
this case, however, the characteristic functions are 

(3.16) P(z,/c) = 1 - z -2 - Xaz _1 (/c - /c' 1 ) 
and 

(3.17) R(z,k) = 1 - E (/) (-1) /+1 * -/ K y = (1 - z~ 1 k Y . 


Therefore, 0(z,ac) > 0 at all points (z,#c) in (3.4), except for 


fi (z = -1, k — -1) = | P(z = -1, k — -1)| + | R(z = -1, ({sr-l)| sO; 
so (3.4) fails, and Theorem 3.1 does not imply stability. 

This counterexample -- which shows that our criteria in Theorems 3.1 and 3.1' are 
sufficient but not necessary for stability — is not unique. As we shall see, there are 
other well-known approximations for which H(z,/c) does not vanish on the (z,/c) 
domain in (3.4), with the exception of the particular point (z,/c) = (-1,-1). With this 
in mind, we shall prove in Section 4 the following alternative to Theorem 3.1, where 
the point (z,k) — (-1,-1) is treated seperately. 
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THEOREM 3.2 (2nd Main Theorem). Approximation (2.8) is stable if 


(3.18) 


BP[z,k) j dP{ z,k) 
dz f die 


<0 y 

2 =K =-l 


and 


(3.19) 


n{z,k) > 0 for all 

{M = |k| = 1, (z,/c) ± (1,1)} U {|*| > 1, 0 < \k\ < 1} . 


Having Theorem 3.2, we follow the idea that led to Theorem 3.1', and split the 
(z, k) domain in (3.19) into three disjoint sets in order to obtain: 


THEOREM 3.2' (Theorem 3.2 restated). Approximation (2.8) is stable t/ (3.18) 


holds and 




(3.20a) 

Q(z,k) > 0 

for all 

\z\ = |/c| = 1, k 1, (z,k) (-1,-1) , 

(3.20b) 

Cl(z,K=l) > 0 

for all 

\ z \ = 1, 2 7^ 1 , 

(3.20c) 

n(2,/c) > o 

for all 

\z\ > 1, 0 < |/c| < 1 . 


Since the inequalities in (3.20b, c) coincide with those in (3.5b, c), then Lemma 3.1 
immediately yields: 


LEMMA 3.2 (2nd Main Lemma). Inequalities (3.20b) and (3.20c) hold under the 
hypotheses in parts (ii) and (iii) of Lemma 3.1, respectively. 

We deliberately choose not to apply part (i) of Lemma 3.1 to the inequality in 
(3.20a), since if the basic scheme or the boundary conditions are dissipative, then by 
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Assumption 3.3 (i) and Definition 3.2, we have either (3.11) or 


R{z,k) 0, \z\ > 1, \k\ = 1, k 1 . 

So 

D(2 = -1, k = -1) > 0 ; 

and Theorems 3.2 and 3.2' lose any possible advantage they may have had over 
Theorems 3.1 and 3.1'. 

EXAMPLE 3.6. Consider any basic scheme of form (2.8a) whose characteristic 
function satisfies (3.18) as well as 

(3.21) P(2 ,k) 7^ 0, z = k, \k\ = 1, k ± 1 ; 

and determine the boundary conditions by oblique extrapolation as in (3.7). The boun- 
dary characteristic function, given in (3.17), satisfies 

(3.22) Q(z,k) > |P(2 ,ac)| >0, z k . 

In addition, by (3.21), 

(3.23) n(z,/c) > |P(2,/c)| >0, z = k, |k| = 1, k ± 1 . 

So (3.22) and (3.23) imply (3.19), and Theorem 3.2 yields stability. 

The inequalities in (3.18) and (3.21) are met by a number of well known schemes, 
including: 

(a) The unitary Leap-Frog scheme in (3.15) (compare [5, Theorem 6.2)), whose 
characteristic function in (3.16) satisfies 
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dP(z,K 

dz 


1 / *PM 

/ 3/C 


= ~r~ < ° 

z=«=-i \a 


and 


P{z,k) = (1 - Xa)(l -k 2 ) =£ 0 , z = k, k ^ ± 1 , (0 < Xa < 1) 


(b) The unitary, five-point Leap-Frog scheme (e.g. [8,9]): 


(3.24) 


Here, 


»*(' + A/) = «v(f - A() + ^ Ivh(') - v„-i(OI - Y [»„«(') - K-M \ , 


0 < \a < — , 1 / = 2,3,* 
5 


P(z,k) = 1 - z 2 + -g- Xaz _1 (/c - /c ! )(ac + /c 1 - 8) ; 


so 


and 



.P(£,/c) = (1 - K 2 )[1 + — \a(K + K~ l - 8)] 7^ 0 , z = k, |/c| = 1 , K ± 1 . 

6 

(c) The almost-dissipative Lax-Friedrichs scheme (e.g. [10]): 
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(3.25) 


7 


v u {t + M) = ~(vm(<) + + y- k+i(0 - vi(OJ 

0 < \a < 1, v = 1,2,3,* * ; 

for which 


? -i 

(3.26) P(z,/c) = 1 - - 7~ [« + /c" 1 + Xa(/c - /c -1 )] , 

so again, 


(3.27) 


and 



F(2,/c) = ^ (1 ~ Xfl )(l ~ K 2 ) 7 ^ °> 2 = K i M = 1> «^±1. 


We realize that the stability question for these three schemes, combined with 
oblique extrapolation at the boundary, could not have been handled by Theorem 3.1 
since here, 


fi(z = - 1 , ac = —1) = 0 . 

EXAMPLE 3.7 (compare [4, Example 2.7]). Take any basic scheme that satisfies 
(3.18), and let the boundary conditions be determined by the first-order accurate, 
weighted Euler scheme in (3.10). The boundary characteristic function is 

R(z,k) = 1 - z~ 2 - Xa(2Kz~ l - 1 - z~ 2 ) ; 
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so its roots z(/c) satisfy 


(3.28) z(k = «*) = e* , $(£) = s/(\ a ) 2 + e~‘^[l - (Xa) 2 ] . 

A d T 1 

Denoting p = (Xa) 2 , we obtain 

l^(0| 4 = P 2 + (1 ~ P) 2 + 2p(l - p) cos2f< p 2 + (1 - p) 2 + 2p{l - p) = 1 

0 <|e| <n,(0 < p <1) 

hence by (3.28), 

(3.29) M* = e*)l < j^ 11 <1, 0 <K| O. 

Consequently, 

R(z,k)^ 0, \z\ > 1, |/c| = 1, k 7 ^ ± 1 ; 

so 

fi(z,/c) > |/?(z,/c)| > 0, \z\ >1, |k| = 1, K^± 1 , 

and (3.20b) follows. 

Next, since 

n (z = -1, k = 1) > | R(z = -1 ,k = 1)| = 4Xa > 0 , 
then Lemmas 3.1 (ii) (d) and 3.2 imply (3.20b). 

Further, by (3.29), the roots of R{z,k) satisfy 
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|z(/c = e'*)| <1, |£| < 7T ; 

so our boundary conditions fulfill the von Neumann condition. Moreover, we easily 
verify that 


= 1 -Xuk 7^ 0, 0 < |/c| < 1, (0 < Xa < 1) . 

Hence, Lemmas 3.1 (iii) and 3.2 assure (3.20c), and by Theorem 3.2' we have stability. 

We recall that the Leap-Frog scheme in (3.15), the five-point Leap-Frog scheme in 
(3.24), and the Lax-Friedrichs schemes in (3.25), all satisfy (3.18); so they fit into the 
above example. 

We also point out that the stable approximation (3.25) (3.10) was mistakenly 
declared unstable in [4j. 

EXAMPLE 3.8. Consider the almost-dissipative Lax-Friedrichs scheme in (3.25), 
with the Box-scheme boundary conditions in (3.8). As indicated in (3.27), the charac- 
teristic function in (3.26) fulfills (3.18). Also, the roots of P(z,k ) satisfy 

z{k = e'f) = cos £ + tXasin f ; 
so 

z(k = - 1 ) = -1 

and 

(3.30) \z{k = e *)| 2 = cos 2 e + (Xa) 2 sin 2 e < 1, 0 < |£| < tt, (0 < \a < 1) . 

Thus, 

P(z,k)^ 0, \z\ > l, M = 1, K^l, (z,k) ^ (-1,-1) ; 
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hence 


fi(z,/c) > \P(z,k)\ > 0, \z\ > 1, |ac| = 1, k^I, ( 2 , k) ^ (-1,-1) , 
and we obtain (3.20a). 

In addition, Lemmas 3.1 (ii) (c) and 3.2 immediately give (3.20b). Also, by (3.30) 
we get 

= e'*)| <1, U| < x ; 

so the boundary conditions fulfill the von Neumann condition. And since 

Re r_j(/c) = 1 + Re(/c) + Aa[l -Re(/c)J ^ 0, |k| <1, (0 < \a < 1) , 

then Lemmas 3.1 (iii) and 3.2 yield (3.20c), so Theorem 3.2' implies stability. 
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4. Proof of Main Results. 


In [4] we proved: 

THEOREM 4.1 [4, Theorem 2.1]. Approximation (2.8) is stable if 

n(z,/c) > 0 for | z\ > 0, 0 < |k| < 1, [z,k) yA (1,1) . 

With this we can easily obtain: 

Proof of Theorem 3.1. Since by Assumption 3.2 the basic scheme fulfills the von 
Neumann condition, then the characteristic function in (3.1) satisfies 

P(z,k) yA 0 for \z\ >1, |k| = 1 . 

Thus, we always have 

(4.1) fi(*,/£) > 0, |*| > 1, |«| = 1 , 

which together with the statement of Theorem 4.1, completes the proof. □ 

In order to prove Theorem 3.2 we again consider P(z,k). By Assumption 3.4, for 
each z with |z| > 1, P(z,k) has r + p roots k{z). These continuous roots, which 
play a major role in the stability analysis of Approximation (2.8), have the following 
separation property. 

LEMMA 4.1 ((4, Lemma 3.1]; compare [5, Lemmas 5.1 and 5.2]). For \z\ > 1, 
the characteristic function P{z,k) has precisely r roots k(z) with 0 < |/c(z)| <1, 
p roots with |k(z)| > 1, and no roots with |k(z)| =1. 
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By this lemma, the roots k(z) of P{z,k) split for \z\ >1 into two groups: r 
inner roots satisfying |/£(z)| <1, and p outer roots with |/c(z)| >1. Using a con- 
tinuity argument, we realize that these groups of inner and outer roots remain well- 
defined for \z\ > 1 as well, where milder inequalities, |/c( 2 )| < 1 and |/c( 2 )| > 1, 
hold, respectively. Here, of course, if for z 0 with | 2 0 | = 1 we have |k( 2 0 )| = 1, 
then k(z 0 ) is an inner (outer) root of P(z 0 ,k) if and only if, as z, \z\ > 1, 
approaches z 0 , there exists an inner (outer) root k(z) of P(z,k) that satisfies 
/c(z) -+ /c(z 0 ). Since by Assumption 3.4, k=0 is not a root of P{z,k) for \z\ > 1, 
we summarize the above argument as follows: 

LEMMA 4.2 [4, Lemma 3.2]. For \z\ > 1, the r + p roots /c(r) of the basic 
characteristic function P(z,k) split into r inner roots with 0 < |k(z)| < 1 and p 
outer roots with | k (2)| > 1. 

We quote now the following stability criterion: 

THEOREM 4.2 [4, Corollary 3.1]. Approximation (2.8) is stable if and only if for 
every z with \z\ > 1 and each corresponding inner root k(z) of P(z,k), 

H(z,k(2)) = |P(z,k( 2))( + |i?(2,K(2))| > 0 . 

We shall also need: 

LEMMA 4.3. (i) k = 1 is not an inner root of P(z,k) for z — 1. 

(ii) //(3.18) holds, then k = -1 is not an inner root for z = -1. 
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Proof, (i) The first part of the lemma is precisely the content of Lemma 3.3 in [4j. 
(ii) If P(z = -1, k — -1) 0, then there is nothing to prove, so suppose that 


(4.2) P(z = - 1 , k = -1) =0 . 

By (3.18), 


(4.3) 


8P{ Z,k) 
8k 


' Z ==J K 1 


7 ^ 0 ; 


thus having (4.2) and (4.3), we employ the implicit function theorem and find that in 
some neighborhood N of z = - 1 , the characteristic equation 


P(z,k) = 0 

has a unique differentiable solution k = &(r) satisfying 


(4.4) k(z = -1) = -1 . 

Since 


P[z,k{z)) = 0 for z € N , 


then 


(4.5) 


dP[z,K) , 9 P(z,k) d&jz) 


dz 


8k 


dz 


= 0 


Z—K = -l 


Hence, (3.18) and (4.5) imply, 


(4.6) 


dtzjz) 

dz 


Z = -\ 


_ 9P(z,k) j 8 P(z,k ) 

dz 1 8 k 


1 2 =«=-! 


= 7 > 0 
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By (4.4) and (4.6), therefore, for z = -1 - e with sufficiently small e > 0, a Taylor 
expansion gives 


k(z = -1 -e) = k(z = —1) — e 


dk(z) 

dz 


+ 0 { £ 2 ) 

Z -- 1 


= -1 - 7C + 0(e 2 ) < -1 . 

Thus, for z in some left real neighborhood of -1 we have 

(4.7) |/c( 2 )| > 1 . 

Since k(z) is a continuous function in N then by Lemma 4.1, tbe inequality in (4.7) 
holds for all z € N, \z\ > 1. That is, k( 2 ) is an outer root of P{z,k) for z € N, 
\z\ > 1. Since in addition, 

— > _i ( 
z — -1 

then by definition, k = -1 is an outer root for z = - 1 , and the proof is com- 
plete. □ 

Proof of Theorem 3.2. By (3.18) and Lemma 4.3, k = 1 and k = -1 are 
excluded as inner roots of P(z,k) for z = 1 and z = - 1 , respectively. Thus, in 
view of Theorem 4.2 and Lemma 4.2, approximation (2.8) is stable if 

(4.8) Q(z,k) > 0 for all \z\ > 1, 0 < |k| < 1, (*,«) 7^ ± (1,1) . 

Since the basic scheme fulfills the von Neumann condition then, as shown in the 
proof of Lemma 4.1, we automatically have (4.1). But (4.1) combined with the 
hypothesis in (3.19) implies (4.8), and stability follows. □ 
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